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Abstract 

In HEP experiments the description of the trajectory of a charged particle is obtained from a fit to measurements in 
tracking detectors. The parametrization of the trajectory has to account for bending in the magnetic field, energy loss 
and multiple scattering in the detector material. General broken lines implement a track model with proper description 
of multiple scattering leading to linear equations with a special structure of the corresponding matrix allowing for a 
fast solution with the computing time depending linearly on the number of measurements. The calculation of the 
full covariance matrix along the trajectory enables the application to track based alignment and calibration of large 
detectors with global methods. 
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1. Introduction 

The trajectory of a charged particle in a homogeneous 
magnetic field neglecting the interactions with the de- 
tector material is described by a helix. In a global co- 
ordinate system (x,y,z) with the magnetic field in z- 
direction it can be parametrized by the inverse momen- 
tum (times charge) q/ p, an angle <po at and the distance 
do to the point of closest approach in the (x,y)-plane, 
the dip angle A to that plane and the offset zo at the point 
of closest approach. Energy loss in the detector mate- 
rial due to ionization or radiation (for electrons) leads 
to a reduction of the momentum. Multiple scattering, 
mainly due to Coulomb interaction with the electrons 
in the atoms, results in random changes in direction and 
spatial position with expectation values of zero and vari- 
ance depending on the traversed material and the parti- 
cle momentum. In addition the magnetic field is usu- 
ally homogeneous only in an approximation. Therefore 
more advanced track models are required. 

The effect of multiple scattering can be taken into ac- 
count in different ways (2. For global methods it can 
be added to the measurement errors leading to a non di- 
agonal covariance matrix or explicitly fitted by scatter- 
ing angles as additional track parameters. In both cases 
this requires the inversion of a large matrix of order n 
(number of measurements or number of scatterers) with 
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computing time 0(n y ). Progressive methods |2] like the 
Kalman filter [ 3 1 updating the track parameters for each 
additional measurement and scatterer with a computing 
time 0(n) have become the standard method for track 
fitting. 

The broken lines method is a fast global track refit 
adding the description of multiple scattering to an ini- 
tial trajectory and able to determine the complete co- 
variance matrix of all track parameters. This allows 
the usage as track model for track based alignment and 
calibration with the global Millepede-II |4| algorithm. 
Corrections and covariance matrices for the local track 
parameters (at single points) can be determined with a 
computing time 0{ri) exploiting the sparsity of the ma- 
trix of the corresponding linear equations system. 

The original formulation [5| describes the case of 
a tracking detector with a solenoidal magnetic field 
and independent two-dimensional tracking in the bend- 
ing and perpendicular to the bending plane. It con- 
structs the planar trajectories from the measurements 
including the material around those as thick scatterers. 
In the presence of measurements with components in 
both planes a single trajectory in three dimensions is 
required. General broken lines are describing trajec- 
tories in space in an arbitrary magnetic field for ar- 
bitrary measurements and material distribution. The 
three-dimensional broken lines trajectories from one- 
or two- dimensional independent measurements dedi- 
cated to the track-based alignment of the CMS silicon 
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tracker [6 | with Millepede-II J4] are an intermediate 
step. They are constructed from the material between 
the measurement planes by equivalent thin scatterers. 
Measurements are parametrized by the interpolation of 
(small) offsets perpendicular to the track direction de- 
fined at those scatterers and a common correction Aq/p 
of the inverse momentum. Multiple scattering kinks are 
defined by triplets of thin scatterers. The corresponding 
linear equation system Ax = b with a bordered band 
matrix A with a size n depending on the number of scat- 
terers n scat , band width m (usually 5) and border size 
b — 1 is built from the Jacobians of the propagation 
of the offsets between measurements and scatterers on 
the initial trajectory and solved by root-free Cholesky 
decomposition. In the following the further features of 
general broken lines are presented, the procedure is dis- 
cussed and a comparison with the Kalman filter is per- 
formed. 

2. Towards general broken lines 

2.1. Local coordinate system 

At each measurement plane and each thin scatterer a 
local (orthonormal) coordinate system (m,U2,w) is de- 
fined. The natural choice for the w-direction is perpen- 
dicular to the measurement plane for a measurement 
and parallel to the track direction for a scatterer. A 
local system moving with the track is the curvilinear 
frame (x_ L ,y ± ,z±) with z± in track direction and x± in 
the global (x, y)-plane. At each thin scatterer the two- 
dimensional offset u = (u\,ui) is defined as fit parame- 
ter. 

The variance V* of the multiple scattering kinks k 
in the local system can be calculated from the variance 
6>y ( q j ) [ 7 ] in the curvilinear system by parameter trans- 
formation |8| from the curvilinear to the local slopes 
which depends on the projections c, = e tlac k • e Uj of the 
offset directions e Ul of the local frame onto the track di- 
rection eu-ack: 
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With at least one offset defined perpendicular to the 
track direction this is a diagonal matrix. 

2.2. Local track parameters 

In the local frame a track can be described by an off- 
set {u\,U2), a slope (u' x ,u' 2 ) = fl( ^" 2) and the inverse 



momentum q/p. The general broken lines fit deter- 
mines from the fit parameters x = (Aq/ p, Ui, .., u, w ) at 
each thin scatterer i track parameter corrections Api oc = 
(Aq/p, u', Uj). The slope before the scattering is defined 
by (Aq/p,Ui-i,VLi) and afterwards by (Ag/p,u,, u,+i) 
(eqn (6) in |6|). At the measurement planes the track 
parameter corrections Api oc = (Aq/p,vL t ,VLi n t) can be 
obtained from interpolation of the enclosing scatterers 
(eqn (8) and inserting U; nt for Uq in eqn (6) in J6)). In 
any case Api oc depends only on Aq/p and two adjacent 
offsets. Therefore for the covariance matrices of the 
track parameter corrections the bordered band part of 
the covariance matrix A -1 of the fit parameters is suf- 
ficient. For a sparse matrix the elements of the inverse 
matrix inside the sparsity pattern can by calculated by 
special methods without those outside that pattern ||9). 
This allows to obtain the bordered band part of A -1 with 
0(n(m + b) 2 ) operations linear in the number of scatter- 
ers. 

2.3. Measurements 

The measurements m are (the residuals with respect 
to the initial trajectory of) arbitrary observables with 
predictions defined by the local track parameters in the 
measurement plane. This can be one or two dimen- 
sional position measurements or track segments con- 
taining slope information. The linearized prediction 
is m(x) = H m x. In case the corresponding covari- 
ance matrices V,„ are not diagonal they should be di- 
agonalized and the measurements and predictions mod- 
ified accordingly. This allows the usage of univariate 
M-estimators for outlier down-weighting and the inter- 
facing to Millepede-II [4 1 expecting independent scalar 
measurements. 

2.4. Iterations 

In case the resulting fit parameters are not small cor- 
rections as assumed for the linearization of the track 
model (eqn (3) in J6 |) the trajectory has to be iterated. 
This requires initial non-zero values ko in the prediction 
of the multiple scattering kinks: k(x) = H^x + ko. The 
expectation value remains (k(x)) = 0. 

2.5. Seeding 

The general broken lines are seeded by an initial tra- 
jectory. Alternatively to the result of a fit of the mea- 
surements (internal seeding) this can be a prediction 
from another part of the detector (external seeding). The 
seeding track parameters at some reference point are 
used for the propagation along the trajectory according 
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to the magnetic field (and average energy loss) to calcu- 
late residuals for the measurements and parameter trans- 
formation matrices. If the track fit has to be iterated to 
account for nonlinear effects the initial trajectory could 
be based on general broken lines itself. On average no 
change of the track parameters at the reference point is 
expected from refitting: (Ap see d) = (H,x) = 0. 

2.6. Minimization 

The fit parameters are determined by minimizing: 

^ meti s 

X 2 (x) = ^ ( H '«.' x - m ') T V~J ; (H m>i x - m,) 
(from measurements) 
Y (H kJ x + k 0A ) T Vjj (H kJ x + k 0;! ) (from kinks) 
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+ (H s x) r V" 1 (H s x) (from external seed) (2) 



The minimization leads to a linear equation system 
Ax = b with x of size n - 2n scat + 1 : 

ft-mcas ft scat ~ ' 

(3a) 

b = £ - £ Hf.V.jk,,; (3b) 
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3. Procedure 

The initial trajectory provides in the local coordinate 
system for the measurements the residuals m with (di- 
agonal) covariance matrix V m , for the scatterers the ini- 
tial (usually zero) kinks ko with covariance matrix V/t, 
optionally the external seed and for the parameter trans- 
formations from one point (measurement, scatterer or 
seed) to the next one the Jacobians T,- = <9pi oc ,i+i /3pi oc ,i- 

The linearized track model of the general broken lines 
is defined at each point by the offset part (d, S, J) ± = 
du/d(q/p, u', u)± of the Jacobians for the transformtion 
to the next and the previous scatterer. For the forward 
propagation they are calculated from the products of the 
point-to-point Jacobians T, for the backward case ad- 
ditional partial inversions (for the offset part) are per- 
formed. In the next step the non-zero elements of the 
matrices H for the linearized predictions are calculated 
from the 2x2 matrices and vectors J ± , W± = iS^ 1 and 
d± (see eqn (7) and (8) in (6)). The measurements m, the 
kinks k, the derivatives from the sparse matrices H and 



the measurement and multiple scattering errors from di- 
agonal covariance matrices V are the ingredients for the 
local fit in Millepede-II J4). Afterwards the linear equa- 
tion system Ax = b is constructed from the individual 
contributions (H r V _1 H, H r V 'm or H r V _1 ko) of mea- 
surements, kinks and optionally the external seed. The 
band part A„ of the matrix is decomposed into a diago- 
nal matrix D and a triagonal band matrix L with unit di- 
agonal (A„ = LDL r ) and the offsets x„ are determined 
from D, L and b„ by solving Lz = b„ and L r x„ = D 1 'z 
0. The band part of A u 1 is calculated according to (9). 
Using block matrix algebra the complete solution x and 
the bordered band part of A -1 are obtained. From this 
corrections and covariance matrices for local track pa- 
rameters at any point (measurement, scatterer or seed) 
can be derived. 

For all this steps (propagation, prediction, construc- 
tion, decomposition, solution and inversion) the num- 
ber of operations is linear in the number of points or 
parameters and a sizable contribution to the total com- 
puting time. As all predictions Hx depend at most on 
three consecutive offsets the matrix A„ has a band width 
m = 3 • 2 - 1 = 5. In case the transformation matrices 
J and S are diagonal the two offset components are de- 
coupled and the band width reduces to m = (3 — 1)-2 = 4 
speeding up the band matrix operations (0{n(m + b) 2 )). 
In addition the construction of A is faster due to fewer 
non-zero elements in H. 



4. Comparision with Kalman filter 

4.1. Mathematical equivalence 

An externally seeded general broken lines trajectory 
with one measurement has only one offset defined and 
the slope has to be used directly as fit parameter: x = 
Apioc.i = (Ag//?,u'j,Ui). The solution of the linear 
equation system is: 

x=A- 1 [H^ 1 V;> 1 ] (4a) 

A" 1 ^[V^+H^V^H,,,,^ 1 (4b) 

This is equivalent to the filtering step of a Kalman filter 
in the weighted mean formalism (equation (8b) in [3 1): 



x k = C 4 [(C^)~ 



C* = 



((t^+H'V^H* 



k r k 
-1 



(5a) 
(5b) 



As the initial trajectory has been based on the track pa- 



rameters from the external seed the prediction x k is 
zero, Cft corresponds to A -1 , C k r l to V, and m k to 
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In general only the covariance matrices for the track 
parameters at single points are calculated by the Kalman 
filter. For the global covariance matrix from all track pa- 
rameters of the trajectory required by global alignment 
and calibration methods an extension of the Kalman fil- 
ter is described in IfTUl . 

4.2. Computational difference 

Track fitting with the Kalman filter algorithm Oil is a 
sequential procedure adding measurements and scatter- 
ers (process noise) to the trajectory one at the time. The 
propagation of track parameters and covariance matri- 
ces from point to point in the filtering step involves mul- 
tiplications of 5 x5 matrices with 0(5 3 ) operations. The 
filtering in the weighted mean formalism is dominated 
by the inversion of 5 x 5 matrices (<9(5 3 ) operations, C 
in equation ( f5b] i) for the addition of a single measure- 
ment. Using the gain matrix formalism requires the in- 
version of matrices of the size of the dimension d of the 
measurements (0(d 3 ) operations) and several multipli- 
cations of 5 x d matrices (0(5 2 ■ d), 0(5 ■ d 2 ) operations). 
The smoothing step to get optimal local track parame- 
ters at all and not only the last point involves inversion 
and multiplications of 5 x5 matrices (0(5 3 ) operations). 
For small d the gain matrix formalism should be fastest. 

The general broken lines can add several (n, between 
one and all) measurements with one fit . The effort is 
defined by the construction of the linear equation sys- 
tem from the non-zero derivatives from the propagation 
to the previous and next thin scatterer and its solution 
using 0(n par ) ■ (m + l) 2 operations with in general band 
width m — 5 and n par = 2n scat + 1, depending on the 
number of thin scatterers n scat = 2..2n. 

In a toy setup the performance of track refitting has 
been compared. Tracks in a detector consisting of 10- 
50 planes with two independent measurements (d = 2) 
and thin scatterers coinciding with the measurement 
planes (simple model for a silicon tracker) or one or two 
thin scatterers between adjacent planes have been sim- 
ulated. The magnetic field B is assumed to be constant 
between adjacent points and the initial trajectory uses 
a detailed helix propagator [8 1 or the simplified version 
for weak deflection (limit \B\/p — > 0) as in [6|. Start- 
ing from the initial trajectory corrections and covariance 
matrix of track parameters at one or both track ends have 
been calculated. The results of this test are shown in fig- 
ure [T] For fit results at one track end Kalman filtering is 
compared with a general broken lines fit including one 
transformation from fit to local track parameters. De- 
pending on the number of scatterers per measurement 
the fit with general broken lines is almost as fast as the 
Kalman filter for the detailed propagator (±25%) and 
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Figure 1: Computing time for track refitting using different algo- 
rithms, propagation Jacobians and setups of scatterers normalized to 
the fastest case. The top plots are with the detailed helix propagator, 
the bottom ones with a simplified version, the left ones for filtering 
(one track end) only and the right ones for filtering and smoothing 
(results at both ends). The methods shown are general broken lines 
(GBL) and Kalman filtering in weighted mean (KWM) or gain ma- 
trix (KGM) formalism. The case of the scattering planes coinciding 
with the measurement planes is indicates by circles, one thin scatterer 
between measurements by squares and two thin scatterers between 
measurements by triangles. 

faster for the simplified one (10-80%). This speedup 
is due to the reduced number of non-zero elements in 
H for the matrix construction and reduced band width 
m = 4. For the Kalman filter the gain matrix formal- 
ism is a little faster than the weighted mean formalism. 
For fit results at both track ends Kalman filtering and 
smoothing is compared with a general broken lines fit 
including two transformations from fit to local track pa- 
rameters. In this configuration the general broken lines 
are faster in all cases tested (30-300%). 

5. Summary 

A trajectory based on General Broken Lines (GBL) 
is a track refit to add the description of multiple scat- 
tering to an initial trajectory based on the propagation 
in a magnetic field (and average energy loss). The ini- 
tial trajectory can be the result of a fit of the measure- 
ments (internal seed) or a prediction from another de- 
tector part (external seed). It is constructed from a se- 
quence of (pairs of) thin scatterers describing the mul- 
tiple scattering in the material between adjacent mea- 
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surement planes. Predictions for the measurements are 
obtained by interpolation between the enclosing scatter- 
ers and triplets of adjacent scatterers define kink angles 
with variance according to the material of the central 
scatterer. The required propagation (on the initial trajec- 
tory) from a measurement plane or scatterer to the previ- 
ous and next scatterer is using locally a linearized track 
model. This linearization may necessitate iterations of 
the fitting procedure. The matrix of the resulting linear 
equations system is a bordered band matrix allowing for 
a fast solution by a root-free Cholesky decomposition 
in a time proportional to the number of measurements. 
At all scatterers and measurement planes corrections for 
the local track parameters are determined. The corre- 
sponding covariance matrices and the pulls for the mea- 
surements and kinks require only the calculation of the 
bordered band part of the inverse matrix. 

In a simulated detector, track refitting with GBL is 
about as fast or a little faster than Kalman filtering, 
depending on the multiple scattering setup and on the 
complexity of the propagation Jacobian. In comparison 
with Kalman filtering and smoothing it is up to a factor 
3 faster. 

A GBL fit with external seed and one (additional) 
measurement is equivalent to the filtering step of the 
track fit with a Kalman filter. 

As the track refit can provide the complete covari- 
ance matrix of all track parameters GBL are well 
suited as track model for alignment and calibration with 
Millepede-II g). 
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